Spin and mass imbalance in a mixture of two species of fermionic atoms in a ID optical lattice 
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In this paper, we study the role of both "spin"(species) and mass imbalance in a mixture of two species 
of fermionic atoms with attractive interaction in an one-dimensional optical lattice. Using the bosonization 
approach, quantum phase transitions between a liquid phase and phase separated states are studied under various 
conditions of interaction, spin imbalance, and mass imbalance. We find that, in the phase-separated region, there 
exists two kinds of phase separation and a special quantum phase transition might exist between them in the large 
U limit. On the other hand, the singlet superconducting correlation dominates in the liquid phase. The pairing 
behavior has been also demonstrated that there is oscillating behavior in real space. We find both the spin and 
mass imbalance are in favor of the formation of Fulde-Ferrell-Larkin-Ovchinnikov state. 
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Ultracold atoms trapped in optical lattices have been in- 
creasingly used to simulate the rich physics in strongly cor- 
related condensed matter systems. In particular, the Fulde- 
Ferrell-Larkin-Ovchinnikov (FFLO) state (H 01 in a mag- 
netized superconductor, in which spin imbalance of elec- 
trons leads to Cooper pair formation and condensate in states 
with nonzero momentum (inhomogeneous distribution in real 
space), was proposed four decades ago. However, the ob- 
servation of the FFLO state in solids has been proven to be 
extremely difficult due to the Meissner effect, and was only 
achieved recently in a heavy fermion system y|] . 

Experimentally, two hyperfine states of ultracold fermionic 
atoms play the roles of up and down spins. Their popula- 
tion could be controlled by using radio-frequency field, such 
that the hope of observing the FFLO state in cold atomic sys- 
tems has been renewed recently 0, H . Since the dimen- 
sionality of the cold atomic systems can be easily tuned, and 
indeed cold atoms have already been successfully trapped in 
one-dimensional (ID) waveguide, it seems natural to consider 
these nonhomogeneous pairing behaviors in this low dimen- 
sional systems. 

So far, much attention has been paid to ID spin-polarized 
fermionic systems with attractive interaction by using differ- 
ent methods and techniques, such as Bethe-ansatz [7], density- 
matrix renormalization-group [8], Quantum Monte Carlo J9[], 
and bosonization II 1011 . All of these studies are based on ap- 
proximating the system by the ID Hubbard model in either 
uniformly distributed or harmonic trap cases respectively. In 
uniformly distributed case, the dominant order in the ground 
state is singlet superconducting (SS) |9|], and existence of 
FFLO state l7ll8ll9l [Ioll . which is oscillating in real space and 
peak in momentum space (Non-zero momentum of cooper 
pairs). As polarization increasing gradually, phase separation 
will happen gradually in the Harmonic trap, like FFLO & fully 
paired wings, FFLO & polarized wings and Chandrasekhar- 
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Recently, Taglieber et al flUl have successfully trapped a 
quantum degenerate Fermi-Fermi mixture, i.e. 6 Li and 40 K, 
by evaporating cooled bosonic 87 Rb gas. The advance raises 
a new and interesting question. What is the influence of mass 
difference to the spin polarized fermionic system or does it 
become easier to observe the FFLO state in such a system. In 
Ref. IU2I1 . the influence of mass difference has been discussed 
with TVf = N[. They found that there is a phase transition 
between SS and charge density wave (CDW) in the negative 
U case, and phase separation exists in the positive U case. 
The bosonization study lfl3ll showed that the phase boundary 
between the density-wave phase and phase separation scales 
like U 2 in the weak coupling region. The results are consistent 
with previous numerical studies [ 14] by the exact diagonaliza- 
tion and density-matrix renormalization group technique. 



The main purpose of this paper is to study the role of 
both "spin"(species) and mass imbalance in a mixture of two 
species of fermionic atoms with attractive interaction in op- 
tical lattices. We use the ID asymmetric Hubbard model 
(AHM) as an effective model to describe such a mixture in 
the optical lattice, and then study the phase separation, domi- 
nant order and pairing behavior in the context of different spin 
populations and negative U region using the bosonization ap- 
proach. We find that if the system is partially polarized, there 
can exist two different phase-separated states, and the FFLO 
state might be more stable in the presence of mass imbalance. 



The paper is organized as follows. In section|IIl we obtain 
the bosonized form of the ID asymmetric Hubbard model, 
then simplify it under conditions of away from half filling 
n < 1 and U < 0. In section Hill we study the ground-state 
phase diagram, and the effects of spin imbalance, and mass 
imbalance. In section[V] we discuss dominant order and pair- 
ing behavior in our interested system. Finally, we summarize 
our results in section [Vll 
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II. THE ASYMMETRIC HUBBARD MODEL AND ITS 
BOSONIZED FORM 

Here we consider a mixture of two species of fermionic 
atoms loaded in a ID optical lattice. In experiments, Such 
an optical lattice potential can be written as 



V(x,y,z) = V sm 2 (kx) + V ± [s™ 2 (ky) + sin 2 (fe)] 



(1) 



action between two species of atoms. In the following we set 
t-\ as unit, N lT = £j rijg-, band filling n = (iVf + N[)/L, mass 
ratio t = ti/tf = m\jmi, mass imbalance z-(t-\- fj)/(tf + f|), 
and polarization P = \N<\ - Ni\/(Nj + N±). 

In the standard bosonization method lfl5l[l6ll . the AHM can 
be expressed in terms of canonical Bose fields and their dual 
counterparts as 



where Vo(Vj_) = Vo(v_l)£r in unit of the recoil energy Er = 
H 2 k 2 /2m. If Vj_ » Vo, the hoping process in the yz plane is 
frozen; while along x direction, the hopping integrals depend 
on the mass of atoms. So the system is quasi one-dimensional. 
Without loss of generality, we use "spin" cr =f, 1 to denote the 
type of atoms. For sufficiently low temperatures, the atoms 
will be confined to the lowest Bloch band, then the system 
can be described by the ID AHM, whose Hamiltonian reads 



(2) 



where ct ff (c JO -) are fermion creation (annihilation) operators 
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, and U the on-site inter- where all parameters take the following forms 
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— = Jl + U - . (6) 

K c V 27r[fj sin(fcf-fa) + fj sin(fc/rja)] 



— = Jl U - , (7) 

K s V 27r[f| sin(fcf-j-a) + f j sin(fc/rja)] 



Av = a[f| sin(fe f |-a) - fj_ sin(fe f |a)]. (8) 

Here the bosonic fields <p c and </> s characterize the charge and 
spin degree of freedom,respectively. kf^ and kf± are the Fermi 
wavevectors for up- and down-spin atoms, which are deter- 
mined by number density of each component, and a is the lat- 
tice constant. v e , are the propagation velocities of the charge 
and spin collective modes of the decoupled systems (Av = 0), 
and K c s are the stiffness constants. 

There are two oscillating terms, Umklapp and backward 
terms, in our low energy effective Hamiltonian. If kpf + kpi + 
n/a or kpf + kp^, they will vanish after performing such inte- 
grals. Physically, quasi-momentum conservation laws do not 



hold in the low energy region for these two processes. Even 
if both of them survive, it does not mean that they will con- 
tribute significantly in the long wavelength scale. Accord- 
ing to renormalization-group analysis fll51|, Umklapp term 
cos( V87r0 f ) contributes effectively only for U > case, there 
will be a gap in the charge excitation spectrum. On the other 
hand, backward term cos( VIt^s) contributes effectively only 
for U < case there will be a gap in spin excitation spectrum. 

Therefore, both of them will disappear in our systems, i.e., 
spin imbalance and attractive on-site interaction U < 0. After 
do one loop approximation, the effective Hamiltonian can be 
written as 



H 



eff 



dx 
dx 



(d xt f> c y + K c n 2 c 



+ Av J dx [n c n s + d x <f> c d x <ps\ 



— (d x 4> s ) 2 + K s n 2 

A, 



(9) 



Here we want to emphasize that the coupling constant Av 



3 



depends on both spin polarization and mass difference. Unlike 
the Hubbard model, here N T > N± and iVj < iVj, will be related 
to different cases because up- and down-spin atoms can be 
distinguished from each other due to mass difference. In the 
following, we try to study the consequences of mass imbal- 
ance and spin imbalance for phase separation and dominant 
orders. 



III. PHASE SEPARATION 



The quadratic effective Hamiltonian [Eq. (O] can be diago- 
nalized in terms of two new fields <p ± which are combinations 
of spin and charge degrees of freedom. The corresponding 
velocities have been obtained 



vl + v 2 



+ Av 2 
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As Av — > 0, v + — > max(v c , v s ), v_ — > min(v f , v s ), and here 
v c < v s . As Av increases, v_ decreases until it vanishes at the 
points: 



Av 



Avf = v c v s K c K s 



K C K S 



(11) 
(12) 



At these points, the freezing of the lower bosonic (mixture of 
real spin and charge) mode is accompanied by a divergence 
in the charge and spin response functions. The static charge 
compressibility k diverges at Av = Avi or Av = Av2. It be- 
haves as 



K= K 



Av 



Av 
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Beyond these points, the susceptibilities become negative. 
This behavior of the static response functions together with 
the vanishing of the collective mode velocity indicates that 
the ground state becomes unstable [ 17] and undergoes a first- 
order phase transition [18]. The instability is known as phase 
separation and has been shown to occur in the extended Hub- 
bard Model and in the t — J model lfl9l l20ll . In our case, we 
obtain 



Avi = a(f| sin(A:/r|fl) + f j sin(^a)), 



(14) 



U 2 

Av 2 = + y [a(r T sin(fc FT a) + t i sin(£ n a))] 2 - ^-^-{^) 

It is obvious that if Av 2 > Av?, the system is in PS region. 
After doing some calculations IU3I1 . we arrive at the condition 
of phase separation: 
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It seems that up- and down-spin are symmetric in above ex- 
pression (Bosonization method only work in the weak cou- 
pling region). However, it is not the case in the large U limit. 
In the following, we will discuss these two cases in more de- 
tails. 
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FIG. 1: The phase boundary between the phase separation (below the 
line) and liquid phase (above the line) in the t -\U\ plane. Here the 
polarization P = 0.5, band filling n = 0.8, and U < 0. 



In negative U case, atoms with opposite spin try to form 
singlet pairs and lower the ground-state energy further. Sin- 
glet pairs can be regarded as a kind of quasi-particle of mass 
/«! + m±. Generally speaking, there are three kinds of parti- 
cles, up- and down-spin atoms, and bound pairs. The reason 
for phase separation is the large difference of their mass, i.e., 
heavy atoms will stay together and give more space for other 
atoms to hop more freely. In Fig. [1] we show the bound- 
ary of phase separation in the weak coupling limit based on 
the Bosonization approach for fixed polarization P = 0.5, the 
more large interaction U is, the more pairs are (more unoc- 
cupied sites); and the heavier down-atoms are, the difference 
of mass of these three kinds of particles will large. Both of 
these two situations will lead to phase separation easily. Here 
the phenomena we obtained is very similar to observations in 
experiments El lU • 

The phase boundary in the t — P plane under fixed interac- 
tion is obtained similarly. If the polarization is larger under 
fixed band filling, there will be less down-particles, and the 
system will be easily in phase separation region composed of 
up-particles and bound pairs. The phase transition boundary 
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that there are two possible configurations (See Fig. |4|. 
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FIG. 2: The phase boundary between the phase separation (below the 
line) and liquid phase (above the line) in the t — P plane. Here band 
filling n = 0.8. 
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FIG. 3: The phase boundary (solid line) between the phase separation 
(top-right region above the solid line) and a liquid phase (bottom-left 
region below the solid line) in the P — \ U\ plane for t = 0.5 (a) and 
t = 0.2 (b). The Bosonization method is only valid on the left side of 
the dashed lines. Here the band filling n = 0.8. 



is shown in Fig. [2] Two phase diagrams in the U - P plane 
for both cases of t — 0.5 and t = 0.2 are shown in Fig. [3] In 
both cases, if the polarization P is large, the phase separation 
state becomes more stable. If we compare two cases, the mass 
imbalance will also strongly affect the phase boundary, which 
will dragged to the left as the mass imbalance becomes larger. 



IV. PHASE SEPARATION IN STRONG COUPLING LIMIT 

Bosonization approach can not provide explicit informa- 
tion in the phase separation region except for indicating the 
phase boundary. Therefore, we try to obtain some insights 
about configurations in phase separation region by analyzing 
the strong coupling limit, i.e. U — > — oo. In this case, all of di- 
lute particles are paired with their partner. It is easy to realized 
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FIG. 4: Two possible configurations in strong coupling limit, one is 
pairs staying together and unpaired particles moving; another one is 
dilute particles moving in the background of dense particles that are 
congregated together. 
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FIG. 5: The ground-state energy as a function of the polarization for 
different dominant configurations: (a) more light atoms and (b) more 
heavy atoms. Here band filling n = 0.8, /f = 1, fj. = 0.5. 

It is obvious that there are two cases, more light particles 
and more heavy particles. Now we consider the former firstly. 
In this case, there are two dominant configurations for the 
phase separation, i.e. the configuration I of all atom pairs stay- 
ing together and unpaired light atoms moving in free space, 
and the configuration II of heavy atoms staying together and 
light atoms moving in the background of heavy atoms (See 
Fig. |4|. Since the interaction is infinite (only kinetic energy- 
just tight bonding model), the ground-state energy of the two 
configurations can be calculated by 



Er - -r 



1 - 



cos ^TT - COS N-N l+ l 



1 - COS 



cos-^-cos Wt+1 



1 - COS 



AT T + 1 



respectively. We show the ground-state energies for both con- 
figuration, i.e. E(l) and E(IT), as a function of the polariza- 
tion in Fig. |3a). From the figure, we can see that E; is al- 
ways smaller than En- Therefore, the configuration I is the 
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ground state. It can be understood in the following way. For 
light atoms, the hoping integral is large, so they can effectively 
lower the ground-state energy. While, in the configuration II, 
motion of heavy particles contribute less to the ground-state 
energy. So the configuration I is favorable. 
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FIG. 6: Phase diagram for band filling n = 0.8 . It is configuration II 
below the curve and configuration I above the curve. 

On the other hand, if the heavy particles is dense, then 
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The two energies are shown in Fig. EJb). In this case, if 
P < 0.65, the ground state is dominated by configuration II, 
while if P > 0.65, it is by configuration I. In configuration I, 
more heavy particles will contribute to the total energy, even if 
the value of contribution of each heavy particle is smaller. In 
the configuration II, value of contribution to total energy from 
each light particle is larger contribute to the total energy, but 
the number of light particles is smaller. Therefore, there is a 
transition point between two phase separated states. We show 
the phase diagram in Fig. [6] Here we would like to point out 
that the two phase separated states, i.e. I and II, are qualita- 
tively different. For the configuration I, the pair-pair correla- 
tion function (n,|n,|«, + ,i-n, + 4) has a non-vanishing long-range 
behavior. So it is a true long-range order. While for the con- 
figuration II, the correlation function decays algebraically. 



of spin and mass imbalance to the pattern in ground state. For 
this purpose, we need to calculate the correlation functions 
related such kinds of order, like CDW, SDW, and SS. They 
are defined as lfl5ll 



where 



R v (x, x') = (O v (x)0'l(x') + h.c.) 
Ocdw(x) = J^e-^^x^ix), 



Oss(x) = ^ 



i(k Ftr -k F -„.)x , t 



(17) 

(18) 
(19) 
(20) 



for CDW, SDW, and SS, respectively. In ID Fermi liquid, it 
is well known that there is no true long-range order, and the 
correlation function usually has the behavior of a power-law 
decay. 



R v (x, x') oc 



1 



(21) 



In Fig. |7] we show the dependence of three correlation ex- 
ponents on P, \U\, and z respectively. From the figure, we can 
see that the SS correlation dominates in all three cases. Unlike 
the case in Ref. [12], there is no phase transition from SS to 
CDW now. We believe that the system will come into phase 
separation region before this transition happens. 

At the same time, we can obtain oscillating behavior of 
pairs very easily. In our system, pair correlation mainly comes 
from singlet pair correlation because of on-site attractive in- 
teraction. Therefore, SS correlation function is just expected 
now. It is obvious that pair correlation function will take the 
following form: 



Rss(x) = 2 cos [(fc FT - k Fl )x] 



1 



(22) 



If and only if + Ni, the pair correlation will take the os- 
cillation form, and different mass will modify the correlation 
exponent ass together with spin imbalance. This state has 
non-homogeneous distribution in real space, and it is called 
FFLO state. From Fig. 13 we can see that the exponent ass is 
suppressed, hence the FFLO state becomes favorable, as the 
mass difference increases. Therefore, we believe that the mass 
imbalance is good for the FFLO state. 



VI. SUMMARY 



V. DOMINANT ORDER AND PAIRING BEHAVIOR IN 
THE LIQUID PHASE 

As we known, the ground state is SS in standard Hubbard 
model with negative U. In Ref. there is phase transition 
from SS to CDW driven by different mass in negative U case 
but with Ni - N[. And now we want to find the modification 



In this paper, from the bosonized form of the ID asymmet- 
ric Hubbard model, we have studied the role of spin imbal- 
ance and mass imbalance at zero temperature under various 
conditions. The conditions of phase separation have been pre- 
sented, and the effects of spin imbalance, different mass and 
interaction are also discussed in detail, we find that more- 
light-particle case is same as the more-heavy-particle one at 
the weak coupling limit. And they are different in the strong 
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FIG. 7: Three exponents as a function of P (left), \U\ (middle), and z (right). Here band filling n = 0.8, /j = 0.5, V = -2 (for the left one and 
right one). The SS order is always dominating in the liquid phase. 



coupling limit. In the liquid phase, SS is always dominant and 
there is no phase transition from SS to CDW. Finally, our re- 
sults show that the mass imbalance might be in favor of the 
FFLO state. 
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